TM  No. 
791218 


NAVAL  UNDERWATER  SYSTEMS  CENTER 
NEW  LONDON  LABORATORY 
NEW  LONDON,  CONNECTICUT  06320 


Technical  Memorandum 


APPLICATION  OF  LINEAR  PREDICTIVE  SPECTRAL  ANALYSIS 
TO  MULTIPLE  TONES  IN  NOISE 


DATE:  12  December  1979 


*****  dd^r/f/M 

ALBERT  H.  NUTf ALL 
Special  Projects  Dept. 
Advanced  Systems 
Technology  Division 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 


20030930  062 


TM  No. 
791218 


ABSTRACT 

An  application  of  linear  predictive  spectral  analysis,  in  particular 
the  forward- and- backward  averaging  procedure,  to  tones  in  noise  is  made.  The 
accurate  tone- frequency  estimation  capability  and  lack  of  tone-splitting  make 
this  an  attractive  technique  for  spectral  analysis  and  estimation  of  tone 
amplitudes.  The  effect  of  inadequate  pole-order  on  the  spectral  estimate  in 
white  and  colored  noise  is  investigated  in  a  number  of  examples.  The  examples 
here  indicate  that  the  linear  predictive  technique  should  receive  wider  appli¬ 
cation  for  spectral  analysis  in  the  realm  of  underwater  acoustics .  Two  programs 
in  BASIC  for  the  Hewlett  Packard  9845  for  the  cases  of  real  and  complex  data 
have  been  written  and  are  available  from  the  author;  they  are  modifications  to 
the  original  program  in  FORTRAN  written  by  S.L.  Marple. 
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INTRODUCTION 

The  linear  predictive  technique  of  spectral  analysis  has  received  a  great 
deal  of  attention  lately;  see  refs.  1  and  2  and  the  references  listed  therein.  A 
comparison  of  several  linear  predictive  techniques  was  made  in  ref.  3,  with  the 
conclusion  that  the  then  newly-proposed  forward- and-backward  (F&B)  averaging 
technique  and  the  Burg  technique  (ref.  4)  were  as  good  as  any  currently  available. 
The  reason  that  the  Burg  technique  has  been  preferred  over  the  F&B  technique  is 
not  due  to  the  quality  of  the  spectral  estimate  (see  ref.  3,  page  81),  but  main¬ 
ly  that  a  recursive  procedure  exists  for  the  solution  of  the  matrix  equations 
for  the  Burg  technique,  whereas  a  full-order  matrix  inversion  was  required  for 
the  F&B  technique  at  each  filter-order  considered.  The  fact  that  the  F&B  technique 
can  occasionally  (rarely)  yield  an  instable  correlation  recursion  or  a  spiky 
spectral  estimate  (ref.  3,  pages  41  and  35) ,  is  not  a  severe  limitation  of  the 
technique  in  practice;  in  fact,  for  data  sets  of  100  points  or  more,  this  author 
has  not  yet  seen  any  examples  of  poles  outside  the  unit  circle  in  the  z-plane, 
although  some  have  come  extremely  close. 

A  recent  contribution  by  Marple,  ref.  5,  has  now  altered  the  computational 
situation  considerably.  Specifically,  he  has  derived  a  recursive  algorithm  for 
the  solution  of  the  matrix  equations  of  the  F&B  technique ,  where  the  number  of 
computations  are  comparable  to  that  required  for  the  Burg  technique.  Furthermore, 
he  has  shown  that  the  F&B  technique  is  not  subject  to  line-splitting  (ref.  6),  as 
the  Burg  technique  is,  nor  are  the  spectral  peaks  as  biased  when  estimating  the  tone 
frequency  location.  The  reason  for  these  desirable  features  is  that  the  F&B  tech¬ 
nique  is  an  unconstrained  p-th  order  least  squares  minimization,  whereas  the  Burg 
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technique  is  a  constrained  progressive-minimization  procedure,  structured  so  as  to 
satisfy  a  preconceived  recursive  algorithm  (ref.  3,  sections  4.6  and  4.7). 

Accordingly,  it  is  now  reasonable  and  well-advised  to  reconsider  the  FSB 
linear  predictive  technique  for  spectral  analysis  of  typical  underwater  acoustic 
waveforms  which  can  include  several  tones  of  widely  different  levels  in  colored 
noise.  In  view  of  the  "global"  optimization  property  (ref.  1,  page  572)  of  the 
linear  predictive  technique,  such  high  dynamic-range  spectra  should  not  be  too 
adverse  cases  for  consideration.  This  memorandum  will  consider  several  such 
examples  and  point  out  some  important  features  and  properties  of  the  linear  pre¬ 
dictive  FSB  technique. 
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POWER  CONTRIBUTION  FROM  ONE  POLE-PAIR 
The  power  density  spectral  estimate  of  the  linear  predictive  technique 
takes  the  form,  at  p-th  order,  of  (ref.  1,  page  568) 


?(F)« 


A  E 


A  E 


J[(])!e.xj.(-i'2irPU)|  j  A  (e*p  (i  2rt*))  | 


(0 


where  polynomial 


5  % 


) 


(2) 


A  is  the  time-sampling  increment,  E  is  the  prediction-error  power  at  p-th  order, 
e^*i,  ««d  ta,jf  are  the  p-th  order  prediction-filter  coefficients.  Generalizing 
(1)  to  the  complex  z-plane,  and  allowing  for  complex  data  and  filters,  we  express 
(1)  as 


a*  ( 'A**) 


(3) 


(4) 


meaning  that  the  spectral  estimate  in  (1)  is  given  by 

(P^exj>(  (2-n-f^)  —  T5  (fV 

Now  if  cm  has  a  zero  at  z^,  which  is  near  the  unit  circle  in  the  z-plane, 
then  has  a  pole-pair  at  the  points  zq,  l^z  *,  which  are  both  close  to  the  unit 

circle  and  located  on  the  same  radius  vector;  see  figure  1.  This  gives  rise  to  a 


Q  =  (t 
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peak  in  the  spectral  function  P(f),near  the  frequency 

t  artfe)  amfe) 

*0  “  2 A  IT  'N  IT  7  (5) 

where  4  -  WT  is  the  Nyquist  frequency.  And  if  W*i,  the  spectral  peak 


very  large  and  narrow,  with  3  dB  bandwidth  B  approximately  equal  to 

_  _  2  0-1*0  j_  2  (I  - 1*1)  f 

i  V  “  TP  24  TT*  w  (6) 

If  the  values  of  f,  at  which  P(f)  is  calculated,  have  an  increment  larger 
than  bandwidth  (6),  the  peak  of  this  spectral  component  can  be  missed,  and  a  mis¬ 
leading  spectral  calculation  can  result.  One  way  to  alleviate  this  situation  is 
to  calculate  the  power  contribution  due  to  the  pole-pair  in  figure  1,  and  output 
this  information  in  addition  to  calculated  values  of  the  spectral  estimate.  We 


evaluate  this  power  as  follows:  express  (3)  as 


f  (*)  - 


<R(0 


*  £W 


for  z  near  z  ;  see  figure  1.  The  approximation  to  the  spectrum  in  the  neighborhood 


of  Cq  in  figure  1  is 


ft  (e*pM)  = 


The  area  under  this  approximation  is  called  the  power  contribution  of  the  pole-pair 
corresponding  to  i0)  and  is  given  by  the  integral  over  the  entire  frequency  interval 

■  ....  X  JtSL,  <9) 

where  §  denotes  counter-clockwise  integration  around  the  unit  circle  in  the  z-plane. 
Result  (9)  presumes  |^J<  I  ;  if  |^|  >  |  ,  we  get  jj. 
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Now  we  eliminate  unknown  <K(q)  from  ( 9)j  by  observing  from  (7)  that 


p/r\ _ 


Then  (9),  (3),  (2),  and  (5)  yield  for  the  pole-power 


■p  =  £isl  j±d  = _ E- 


i-i^i 


*  i+W  "  |2AQt  l+l^ 


(10) 


A  e*j> 


(-i  2ir-"£  A  l<)J 


\  -tel 


(11) 


This  calculation  requires  that  one  evaluate  both  the  magnitude  and  argument  of  the 
pole  location  zq  of  #>(*)  under  investigation,  since 


G 


e* 


|>  (i  2-rr-fe4  =  exj>  0  W). 


(12) 


Recall  that  fa]  are  the  filter  coefficients,  and  E  is  the  prediction-error  power, 
all  quantities  evaluated  for  p-th  order.  Alternatively,  z^  is  a  zero  location  of 
OM  defined  in  (2). 

For  a  real  process ,  pole-pairs  of  m  occur  in  conjugate  symmetry,  and  the 
spectrum  P(f)>  given  by  (4)}  is  even  about  f=o.  It  is  then  convenient  to  plot  only 
positive  frequencies,  from  0  to  f  N ,  and  double  the  results  above.  Then  we  have  for 
the  pole-power  contribution 

T0  (/enl  <Jo4)  =  J—  ~  » u  |  J  '  (13) 

'  l^ex^UTri^)!  I+I*1 

The  3-dB  bandwidth  is  unchanged  from  (6) ,  however. 
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VOLTAGE  AMPLITUDE  OF  A  POLE-PAIR 

Although  the  pole  locations  of  ffc)  in  (3)  give  a  good  indication  of  the 
frequencies  of  pure  tones,  the  estimated  power  contributions  due  to  these  poles 
are  not  always  good  indicators  of  the  amplitudes  of  pure  tones,  because  a  slight 
radial  movement  of  a  pole  near  the  unit  circle  causes  a  large  change  in  the  esti¬ 
mated  power  contribution;  see  (11).  This  becomes  a  very  pronounced  effect  for 
tones  with  a  large  SNR  (signal-to-noise  ratio) ,  a  condition  where  we  should  expect 
better  indicators  rather  than  poorer  ones.  Therefore,  an  alternative  technique  to 
using  the  p«\e  -power  contribution  is  preferred. 

Since  we  have  a  good  estimate  f  of  a  tone  frequency  from  (5) ,  we  can 

estimate  the  corresponding  tone  amplitude  according  to 

N  | 

JZ*"  e*f(-l  2-rr£AK?)|  (14) 

H5|  •  > 

j  ,N 

where  i**),  is  the  available  data.  (For  real  data,  we  must  take  2  times  this 
result).  As  a  refinement,  we  could  perturb  f  slightly,  to  find  the  local  maxi¬ 
mum  of  this  function,  ffiis  has  been  found  to  be  a  good  technique  in  practice, 
for  high  SNR  tones.  It  gives  a  low- amplitude  indication  for  spectral  peaks  due 
to  spurious  poles,  and  a  high— amplitude  indication  for  strong  tones  that  are 
actually  present.  We  needn't  perform  a  fine-grained  search  over  all  frequencies, 
because  we  simply  choose  frequency  f  according  to  (5).  Of  course,  we  have  to 
solve  for  the  corresponding  zero  zq  of  ;  however,  a  recursive  procedure  is 

described  for  this  zero  location  in  appendix  A. 
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LINEAR  PREDICTION  OF  TONES 

We  will  concentrate  here  on  multiple  tones  sans  noise .  Anticipation  of 
good  approximations,  via  linear  predictive  techniques  applied  to  tones,  can  be 
justified  by  consideration  of  the  following.  Let  complex  waveform  x(t)  be  com¬ 
posed  of  K  complex  exponentials  of  arbitrary  amplitudes ,  frequencies,  and  phases , 
but  no  additive  noise : 

x(t)  =  2EA  (i 2 H*  + « @k).  (15) 


Then  data  value 


X*  s  x(ha)  s  n  +  i^n)  = 

V  a  i 


(16) 


It  follows  that 


x*-,  - 


(17) 


x»-k  =  ^pC-xp(-;2TTfK^  K)  gk(n). 


Now 


these  K  linear  equations  can  be  solved  for  l£ksK»  and  the  results 


substituted  in  (16).  The  end  result  is  an  exact  expression  for  data  value  Xn  in 


terms  of  past  data  values  according  to 

X„  =  , 

*  k=l 


(18) 


where  coefficients  involve  Am)  j  )sk,m<  K.  That  is,  Xn  can  be 

perfectly  linearly-predicted  in  terms  of  its  past  K  data  values,  when  there  is  no 
additive  noise.  Thus  pure  tones  are  good  candidates  for  a  linear  predictive  ap¬ 
proach.  Of  course,  since  Xn_K  can  be  expressed  in  terms  of  Xh_K+|  , ,  from 
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(18) *  it  also  means  that  backward  prediction  is  perfect,  in  the  absence  of  noise. 

As  examples,  we  have  the  following: 
for  K=1 , 

~  ei  •)  s  e»  j 

for  K=2 , 

*  fo+OV.  ~  (^i  > 

*-  =  fe+e»f  >u  j 

for  K*3, 

\  -  (  «)  Ca+  4.  e3ej)  +  (r>e>e>)  x„.,  j 

where 

e*  *  exj>(i2irfKi).  (22) 

The  general  pattern  for  any  K  is  now  obvious  from  (19)- (21). 

It  should  be  pointed  out  that  if  a  dc  component  is  present  in  the  data,  this 
corresponds  to  one  of  the  .  0.  Alternatively,  if  a  process  contains  no  dc 
component,  but  we  subtract  its  (non-zero)  sample  mean,  we  are  creating  a  dc  com¬ 
ponent  in  the  process,  thereby  requiring  one-larger  order  for  perfect  prediction. 
This  behavior  is  illustrated  in  appendix  B. 


(20) 


(21) 
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EXAMPLES 

Our  first  example  is  five  equal-amplitude  complex  exponentials  in  no 
noise.  As  shown  in  (16)- (18),  perfect  prediction  is  possible  with  a  p»5th  order 
filter,  (except  for  calculator  roundoff) .  We  selected  frequencies 

4“  =  -  ,£5",  -3F,  .22,  JO,  .  yy.  (23) 

"V 

The  results  of  the  spectral  analysis  procedure*  for  N=  100  complex  data 
points  and  a  1024-point  FFT  over  the  ±  Nyquist  frequency  band  are  given  in  dB  in 
figure  2  for  p  =  5.  All  the  plots  are  normalized  so  that  their  peak  value  reaches 
the  top  of  the  70dB  range  plotted.  The  magnitudes  and  arguments  of  the  5  zero 
locations  of  d  (?)  are  given  in  table  1. 

Magnitude 
.999714 
.999822 
.999916 
.999105 
.999013 

Table  1.  Zero  Locations  of  m  for  Five  Tones;  p  =  5. 

It  will  be  observed  that  all  zeros  lie  just  inside  the  unit  circle,  rather 
than  on  it,  and  that  the  frequency  estimate  associated  with  each  zero  is  accurate, 
but  not  perfect.  This  is  due  to  the  subtraction  of  the  sample  mean;  see  appendix 
B.  The  spectral  amplitudes  depicted  in  figure  2  differ  amongst  themselves  by 
12  dB,  even  though  the  true  tone  amplitudes  were  equal;  there  are  two  reasons  for 
*The  sample  mean  was  subtracted  from  the  data  set;  see  appendix  B  for  a  discussion 
of  this  effect.  All  of  the  results  to  follow  also  correspond  to  subtraction  of  sample 
mean. 


Argument /tr  = 

-.649911 

-.379931 

.220052 

.699925 

.850057 


13 


TM  No. 

791218 

this.  One  is  that  the  spectral  function  P(f)  in  (1)  was  calculated  at  only  1024 
points  in  the  iNyquist  frequency  ran^s,  and  therefore  the  very  narrow  widths  of 
the  peaks  in  P(f)  have  not  been  adequately  sampled,  in  fact,  the  3-dB  bandwidths, 
B/fNj°f  the  5  peaks  in  figure  2  are,  from  (6),  (1.82,  1.13,  0.54,  5.79,  6.29)10~4 

respectively,  which  are  smaller  than  the  frequency  sampling  increment  2/1024  = 

-4 

20  x  10  .  This  frequency  sampling  computation  effect  will  contribute  greatly  to 

the  apparent  fluctuation  in  peak  spectral  indication,  and  can  be  alleviated  by 
taking  larger-size  FFTs,  or  by  calculating  P(f)  more  finely  just  in  the  neighbor¬ 
hood  of  the  peak  indication. 

The  second  reason  is  that  the  peaks  of  P(f)  do  not  necessarily  have  to 
reflect  true  tone  powers;  rather  it  is  the  area  under  each  peak  which  is  indicative 
of  that  tone  power;  see  ref .  7  .  A  calculation  of  these  pole-powers  according  to 
(11)  yields  values  .971,  .967,  .945,  .999,  and  .994  respectively,  compared  with 
the  exact  value  of  1.  The  tone  amplitudes ,  as  calculated  according  to  (14),  are 
.970,  .980,  .980,  .981, and  .971  respectively,  instead  of  1.  Perturbation  of  f 

o 

in  (14),  about  the  value  (5),  did  not  significantly  change  the  peak  of  (14)  . 

When  the  order  of  the  predictive  filter  is  taken  to  be  p=4,  instead  of  the 
correct  value  of  5  for  this  complex  data  example,  the  results  are  as  indicated  in 
figure  3.  The  two  clc6.es t  tones  at  .7  and  .85  are  not  resolved.  The  pole— powers 
are  given  as  .99,  1.24,  1.03,  and  1.85  respectively,  indicating  a  lumping  together 
of  the  two  closest  tones .  The  tonal  amplitudes  are  grossly  in  error  since  none 
of  the  5  frequencies  were  estimated  accurately  enough.  This  example  points  out 
the  deleterious  effect  of  too  low  an  order,  f,  for  the  predictive  filter. 
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Our  next  example  consists  of  the  same  five  complex  exponentials  with  white 
noise  added,  to  yield  a  signal-to-noise  power  ratio  (SNR)  of  7.8  dB  per  tone, 
as  measured  in  the  Nyquist  frequency  interval  }  ~FU).  First,  in  figure  4  is 
presented  the  spectral  plot  for  a  fifth -order  prediction  filter,  p=5 .  The  fre¬ 
quency  estimates  according  to  (5)  are  in  error  by  | f/fuj  =.004,  on  the  average. 

The  pole-powers  and  tone  amplitudes  are  1.06,  1.10,  1.01,  .93,  1.22  and  .97,  .88, 

.96,  .67,  1.02  respectively. 

When  we  repeat  this  same  example  with  p  changed  to  10,  the  situation  is 
much  improved,  as  figure  5  indicates.  The  frequency  estimates  are  now  in  error 
by  I-f/fJ  =.0007 ,  on  the  average.  The  pole-powers  and  tone  amplitudes  are  1.05,  .83, 
1.00,  -76,  1.30  and  1.00,  1.04,  .99,  .92,  1.02  respectively.  We  notice  a  signi¬ 
ficant  improvement  in  amplitude  estimation,  which  is  due  mainly  to  the  improved 
frequency  estimates.  The  pole-powers  are  more  erratic  than  for  figure  4. 

When  the  SNR  per  tone  is  reduced  further,  to  0  dB  as  measured  in  the  Nyquist 
band,  the  results  in  figures  6  and  7  are  obtained  for  p=10  and  15  respectively. 

The  pole-powers  vary  significantly  from  the  true  values  for  both  examples;  however, 
the  tonal  amplitude  estimates  for  p=15  in  figure  7  are  rather  good:  1.04,  1.11, 

1.01,  .85,  1.08.  The  frequency  estimates  are  better  for  p=15,  being  in  error 
by  If/fJ  =  .00073  on  the  average.  These  results  are  comparable  to  the  higher 
SNR  results  in  figure  5  where  p=10  was  used.  When  p  is  changed  to  25,  results 
substantially  equivalent  to  figure  7  were  obtained,  except  that  there  were  pro¬ 
portionally  more  peaks  and  sharper  tonal  indications;  see  figure  8.  The  large 
variation  in  peak  values  (21dB)  is  not  too  significant  when  we  recall  the  fre¬ 
quency  sampling  and  tone-power  considerations  discussed  earlier. 

_  —  — 

The  SNR  of  each  tone  is  the  ratio  of  each  tone  power  to  the  total  noise  power, 
in  the  Nyquist  band  (~fN  ,fN) . 
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We  now  consider 

a  case  of  N=1000  real  data 

points  with  9  tones  of  unequal 

value  in  white  noise,  j 

Specifically,  the  tonal  amplitudes,  powers,  and  frequencies 

are  listed  in  table  2. 

Spectral  analysis  results  are  presented  in  figures  9-11, 

Amplitude 

Power 

F  r eque  n  cy /Ny qui s  t 

i 

0.5  • 

.11 

5 

12.5 

.19 

2 

2.0 

.32 

8 

32.0 

.41 

3 

4.5 

.53 

6 

18.0 

.62 

1 

0.5 

.69 

7 

24.5 

.77 

5 

12.5 

.89 

Table  2 .  Nine-Tone  Example 


for  p=18,  20,  and  25,  respectively  for  noise  power  .  Due  to  the  realness  of 

the  data,  we  only  need  plot  the  even  spectrum  from  0  to  f„.  Figure  9,  for 
P=18j barely  gives  any  indication  of  the  weak  tone  at  f/fw  =  .69 ,  despite  the  fact 
that  9  real  tones  can  be  predicted  via  an  18  pole  filter  if  noise  were  absent. 
Increasing  p  to  20  or  25  clearly  indicates  all  9  tones.  The  relative  peaks  in 
figures  9-11  are  not  important,  for  reasons  already  discussed.  In  table  3  are  given 
the  estimated  tonal  amplitudes,  powers,  and  frequencies.  For  all  tones,  the  frequency 


P=20 _  _ p=25 


Amplitude 

Power 

Syq3is?Cy/ 

Amplitude 

Power 

0.9944 

0.524 

.109834 

0.9952 

0.493 

.109844 

4.9931 

11.324 

.190012 

4.9933 

15.868 

.190006 

1.9844 

1.983 

.319871 

1.9903 

2.102 

.319906 

8.0100 

33.175 

.410001 

8.0100 

67.474 

.410001 

3.0033 

4.691 

.529956 

3.0034 

4.015 

.529996 

5.9866 

18.917 

.619992 

5.9869 

13.091 

.619997 

0.9352 

0.495 

.690411 

0.9818 

0.469 

.690132 

6.9967 

23.956 

.770020 

6.9961 

18.266 

.770008 

4.9837 

14.663 

.890013 

4.9849 

18.692 

.890005 

Table  3.  Estimated  Parameters  for  Nine  Tones 
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and  amplitude  estimates  are  better  for  p=25  than  for  p=20;  however,  the  pole  power 

estimates  are  significantly  better  for  p=20  than  for  p=25.  This  appears  to  be  due  to 

the  sensitivity  of  the  pole-power  to  the  proximity  of  the  pole  to  the  unit  circle; 

see  (11)  or  (13).  In  fact,  the  best  pole-power  estimates  for  this  example  were 

obtained  at  p=18.  Thus  there  appears  to  be  a  trade-off  between  the  pole-power  and 

amplitude  approaches,  depending  on  the  proximity  of  each  pole  to  the  unit  circle. 

When  white  noise  power  <T^  is  increased  to  1,  the  value  of  p  must  be  increased 

to  maintain  reliable  results;  an  example  is  presented  in  figure  12.  And  when  the 

noise  power  is  increased  to  10,  the  spectral  estimate  for  p=50  is  given  in  figure  13. 

The  weaker  tones  at  f/f  =  .11  and  .691,  which  have  SNRs  =  -13  dB,  have  erratic  esti- 

N 

mates  for  their  pole-powers  and  amplitudes.  Larger  values  of  p  could  conceivably 
be  advantageous  for  these  lower  SNR  tones .  Since  the  number  of  available  data  points 
is  N=1000  here,  larger  values  of  p  are  certainly  admissable;  in  fact,  Akaike  is 

is 

quoted  (ref.  1,  page  575,  footnote)  as  allowing  p  as  large  as  3N  ,  which  would  be 

about  100  for  this  example  of  N=1000.  The  stronger  tones  in  figures  12  and  13  still 

yield  reliable  frequency  and  amplitude  estimates.  For  example,  in  figure  13,  1  dB 

tone  at  f/f  =  .19  is  estimated  to  be  at  .1901,  with  amplitude  4.91  (rather  than  5) , 

N 

and  the  pole-power  estimate  is  11.2  (rather  than  12.5).  Similarly,  the  equal- 
strength  tone  at  f/f  *  .89  is  estimated  to  be  at  .8902,  with  amplitude  4.69.  The 
correspon dings  estimates  for  the  larger  SNR  case  in  figure  12  were  even  better. 

We  now  consider  some  examples  of  tones  in  colored  noise.  In  particular,  the 
additive-noise  spectrum  is  depicted  in  figure  14  and  is  observed  to  have  a  45"  dB 
dynamic  range  over  the  zero-to-Nyquist  frequency  interval  J  see  appendix  C  for  this  noise 
example.  We  now  add  three  tones  at  the  frequencies 

f  *  .15-4,  .isi,  .Wf*  (24) 
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with  SNRs  -10,  -20,  -30  dB  respectively,  where  these  signal-to-noise  ratios  are 
measured  with  respect  to  the  noise  in  the  entire  zero-to-Nyquist  band.  The  strongest 
tone  at  f  is  located  in  the  region  of  highest  noise  density,  in  fact,  although  the 
signal  strength  drops  by  10  dB  in  going  from  frequency  f ^  to  f ^ ,  the  noise  density 
drops  by  21  dB;  thus  the  local  SNR  fcatio  of  spectral  densities)  is  better  by  11  dB 
at  f2  relative  to  f^  In  going  from  f2  to  f^,  the  signal  strength  drops  by  another 
10  dB,  while  the  noise  density  drops  by  7  dB;  thus  the  " local”  SNR  at  f  is  3  dB 
better  than  that  at  f  .  A  more  thorough  treatment  is  given  in  appendix  C. 

The  spectral  estimate  for  N=1000  and  p=25  is  shown  in  figure  15.  As  anti¬ 

cipated  by  the  "local"  SNR  results  above,  the  sharpest  and  most  accurate  spectral 
peak  is  near  f^.  The  frequency  corresponding  to  the  argument  of  the  pole-pair  is 

written  directly  above  each  peak  in  figure  15  et  seq.,  thereby  enabling  a  ready  measure 

of  accuracy. 

In  figures  16  and  17,  the  data  is  left  unchanged,  and  p  is  increased  to 
50  and  75,  respectively.  The  frequency  estimates  are  progressively  better  in 
figures  15-17.  The  tone  amplitude  estimates  in  figure  17  for  p=75  are  accurate 
within  1.2%,  2.3%,  and  .02%  respectively,  of  the  true  values;  the  pole— powers  axe 
less  accurate,  being  off  by  16%,  16.3%  and  1.7%,  respectively. 

In  the  next  series  of  plots,  the  tones  are  decreased  10  dB  in  strength, 
to  -20,  -30,  -40  dB  SNR  at  f^,.  f ^ ,  f  ,  respectively.  Figures  18-21  correspond  re¬ 
spectively  to  p=75 ,  100,  125,  150. 

The  two  weaker  tones,  of  -30  dB  SNR  at  f/f  =  .45,  and  -40  dB  SNR  at  f/f  = 

N  N 

.65,  are  both  very  well  indicated  for  all  the  values  of  p  considered;  however  the 

strongest  tone,  of  -20  dB  SNR  at  f/f  =  .15,  is  never  really  indicated  at  all.  The 

N 
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only  reason  we  indicated  the  frequency  corresponding  to  the  argument  of  the  pole 

near  f/f  =  .15,  is  that  we  knew  the  true  tone  location  for  this  simulation 
N 

example.  In  practice,  the  spectral  peak  at  f/fN  =  .15  is  not  significantly  larger 
than  its  neighbors  to  be  able  to  declare  it  reliably;  in  fact,  the  spurious  peak 
near  f/f  =  .54  is  larger  than  the  true  tone  at  f/f  =  .15.  This  example  points 
out  that  the  "local11  SNR  (ratio  of  signal  energy  to  noise  spectral  density)  is  more 
important  in  detecting  and  estimating  tones  that  the  "global'*  SNR  (ratio  of  tone  power  to 
noise  power  ).  This  result  is  similar  to  the  standard  FFT  approach.  Also  see  appendix  C. 

The  "ribbon  width”  of  fluctuations,  in  the  spectral  estimates  of  figures 
18-21,  is  uniform  with  frequency,  regardless  of  the  local  density  value;  see  ref.  1, 
page  572.  However,  this  ribbon  width  is  increasing  as  the  order,  p,  of  the  filter 
increases.  This  is  a  manifestation  of  the  resolution- vs.— stability  tradeoff  that 
one  must  accept  in  spectral  analysis.  Specifically,  as  p  increases,  the  frequency 
resolution  improves,  but  the  stability  of  the  estimates  degrades. 

The  relative  amplitudes  of  the  tone  indications  near  f/f^  =  -45  and  .65 
in  figures  18—21  are  not  important,  because  these  peaks  are  narrower  than  the  FFT 
frequency-increment,  and  we  undoubted^ missed  the  maximum  peak  value;  for  example, 
in  figure  21,  the  3  dB  bandwidth  of  the  peak  near  f/f  =  .65  is  .001  whereas  the 
frequency  increment  is  2/1024  *  .002.  Thus,  since  the  pole  shifts  slightly  when  p 
is  changed,  we  can  get  variable  spectral  indications.  The  tone  amplitude  indicators 
are  rather  good  for  figure  21,  being  in  error  by  2.1%,  8.3%,  and  5.6%  respectively. 

The  frequency  accuracy  of  the  two  weaker  tones  is  extremely  good,  despite 
the  low  SNR.  The  minimum  error  is  realized  for  p=125  for  this  particular  example. 
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The  best  value  of  p  to  use  probably  depends  on  N,  the  number  of  data  points  avail¬ 
able,  and  the  local  SNR.  Whether  the  number  of  tones  present  is  a  factor  is  un¬ 
known. 

One  last  example  of  eight  tones  with  various  SNRs  ranging  from  -10  dB 
to  -50  dB  is  presented  in  figures  22  and  23  for  p=150  and  200,  respectively.  The 
frequency  and  SN"R  of  each  tone  is  indicated  under  the  peak  indication,  while 
the  estimated  frequency  is  written  above  the  peak, in  figure  22.  The  two  tones 
corresponding  to  f/f  =  .15  and  .75  are  not  detectable;  on  the  other  hand,  the 
-50  dB  SNR  tone  at  .95  is  clearly  detectable.  Increasing  p  to  200  in  figure  23 
gives  a  rather  clear  indication  of  the  tone  at  .15,  but  the  tone  at  .07  is  marginal 
now;  but  in  neither  case  is  the  tone  at  .75  detected.* 

When  the  strength  of  the  tone  at  f/fN  =  .75  was  increased  by  10  dB,  its 
spectral  indication  increased  by  20  dB.  This  indicates  a  thresholding  effect;  that 
is,  if  the  local  SNR  is  sufficiently  large,  a  good  indication  of  tone  presence  is 
obtained,  but  if  the  SNR  is  decreased,  this  indication  rapidly  disappears* 


*  In  appendix  C,  a  quantitative  measure  of  detectability  is  defined  and  illustrated 
for  this  example  of  eight  tones. 
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The  applicability  of  linear  predictive  spectral  analysis  to  multiple  tones 
in  noise  has  been  demonstrated  in  a  number  of  examples ,  including  some  tones  with 
very  small  SNR.  The  best  filter-order,  p,  to  choose  depends  on  the  number  of 
available  data  points,  N,  and  is  also  probably  a  function  of  the  local  SNR  (ratio 
of  signal  energy  to  noise  density)  at  the  frequency  of  interest.  What  may  be 

indicated  here  is  several  plots,  each  for  a  different  value  of  p,  all  computed  from 
the  same  data,  in  order  to  best  detect  and  estimate  different  level  tones  in  colored 
background  noise.  If  so,  the  recursive  procedure  of  ref.  5  is  very  useful,  since 
one  can  plot  spectral  estimates  at  any  selected  points  of  the  recursion. 

If  only  one  large  value  of  p  is  of  interest,  an  apparent  alternative  pro¬ 
cedure  is  to  use  the  straightforward  matrix  inverse  approach  for  that  order  p  alone. 
This  may  be  quicker  and  more  accurate,  if  the  computational  effort  is  smaller. 

However,  a  pxp  matrix  will  be  necessary  to  store  the  sample  correlations  required 
in  the  FSB  technique,  and  this  can  be  excessive  for  large  p  (>100) .  The  storage 
requirements  of  the  Marple  algorithm  are  very  reasonable,  even  for  very  large  p. 
Nevertheless,  it  would  be  worthwhile  to  conduct  at  least  one  data  set  through  both 
procedures,  in  order  to  ascertain  the  accuracy  and  reliability  of  these  extensive 
data-processing  routines  that  are  being  contemplated. 

Some  results  an  a  comparison  of  linear  predictive  techniques  with  conventional 
Fourier  analysis,  and  resolution  capabilities  of  each,  are  presented  in  refs.  8  and  *). 
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PROGRAMS 

A  Fortran  program  for  the  recursive  solution  of  the  predictive-filter 
coefficients,  via  the  normal  matrix  equations  encountered  in  the  linear  predictive 
F&B  technique,  has  been  presented  in  ref.  5.  This  program  has  been  modified  some¬ 
what  here  and  written  in  BASIC  for  the  Hewlett  Packard  9845  Calculator;  for  the 
two  cases  of  real  data  and  complex  data.  Modifications,  leading  to  minimization 
of  execution  time ,  include :  elimination  of  divisions  in  favor  of  multiplications ; 
use  of  auxiliary  variables  to  minimize  the  number  of  repeated  calculations  and  table 
look-ups;  movement  of  some  variables  outside  loops  to  minimize  look-ups;  eliminate 
some  unnecessary  calculations;  and  correction  of  one  programming  error.  These  two 
programs  are  available  from  the  author;  they  have  been  compared  with  the  direct 
matrix  inverse  procedure,  with  11-decimal  agreement  in  results.  A  simple  example 
and  sample  print-outs  are  given  in  appendix  B. 
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APPENDIX  A.  ZERO  LOCATION  OF  $(?•) 


For  a 
derivative  at 


complex  function  f(z),  suppose  we  have  evaluated  f(z)  and  its 
a  point  z ^  fairly  close  to  a  zero  zq  of  f(z).  Then  we  can  express 


f  (z)  -  f  (z^)  +  f‘ 1  (z^)  (z-z^)  for  z  near  z^  . 


(A-l) 


We  can  then  get  an  estimate  of  the  zero  location  of  f(z)  by  setting  (A-l)  to 
zero;  we  obtain  the  estimate 


z  =  z 
2  1  - 


f  (z1) 


f  (z^ 


(f;  (Z;L)  i-  0). 


(A-2) 


This  procedure  can  be  repeated  by  evaluating  f  and  f*  at  point  ;  the  procedure 
converges  if  f*  (z  )  f  0,  that  is,  if  z  is  a  simple  zero  of  f(z). 

For  our  application,  f(z)  is  given  by  polynomial  (2).  It  is  thus  a  simple 
matter  to  evaluate  the  function  and  its  derivative  at  any  point  z.  The  starting 
point  for  this  procedure  is  obtained  as  follows :  we  evaluate  spectral  estimate 
P(f)  in  Cl)  at  a  large  number  of  equally-spaced  points  in  If|<  2A  ’  ^ 
means  of  an  FFT.  We  then  select  particular  frequencies  where  P(f)  is  very  large 
compared  to  neighboring  frequency  points  and  use  these  as  starting  points  in  re¬ 
cursion  (A-2).  That  is,  a  point  on  the  unit  circle  in  figure  1  very  close  to  Q 
is  used  as  starting  point  z^.  Recursion  leads  to  the  closest  zero  of  dfe)  ,  namely 

z  • 
o 
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APPENDIX  B.  FIVE-TONE  EXAMPLE,  PRINT  OUT,  AND  INTERPRETATION 

The  example  of  interest  here  is  the  one  given  in  (23)  ,  namely  five  equal- 
amplitude  complex  exponentials 

X*  =  Xexp(i2TrfKAn)  ,  )ar?iMs,w>,  (B-1) 

where 

JL  s  2Afk  »  "  35>  21>  *7o>  '  85  •  (B— 2 ) 

L  * 

V* 

The  first  method  of  processing  the  data  is:  subtract  the  sample  mean  and 

then  employ  the  F&B  procedure.  The  computer  print-out  results  are  presented  in 

table  B-l,  and  are  interpreted  as  follows.  The  sample  mean  of  the  data  is  complex 

and  is  given  by  the  number  -,02-i,003727.  The  sample  variance  of  the  data  is 

4.8796.  The  least-squares  algorithm  has  several  different  ways  of  terminating  its 

recursive  procedure;  the  criterion  encountered  here  is  status  1/  which  indicates 

that  the  maximum  specified  order  of  5  was  realized  without  satisfying  any  of  the 

IS  .00337.  TVie  best  fi  lter 

stopping  criteria.  The  prediction— error  power  of  the  best  prediction-filteiy^  weights 
are  complex  and  given  by  the  five  pairs  of  numbers  listed  together.  The  variance 
from  the  graph  is  realized  by  numerically  integrating  under  the  spectral  estimate 
(1)  by  means  of  the  Trapezoidal  rule,  at  the  FFT  points  computed,  and  using  the 
prediction-error  power  given  above.  The  locations  of  the  poles  are  specified  by 
their  real,  imaginary,  magnitude,  and  argument/tr  values,  respectively.  The  power 
content  of  each  pole  and  its  3dB  bandwidth  are  evaluated  according  to  (11)  and 
(6).  And  the  amplitude  of  the  corresponding  tone  is  accomplished  according  to 
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(14).  These  are  the  results  quoted  in  table  1  in  the  main  body  of  this  memorandum. 

The  corresponding  spectral  plot  in  figure  2  is  reproduced  here  in  figure  B-l  for 
easy  reference  and  comparison  to  the  following  results. 

It  was  noted  that  subtraction  of  the  sample  mean  from  a  set  of  data  can 
create  an  effective  dc  component  in  the  original  process,  and  thereby  require  an 
additional  order  in  the  filter  for  its  prediction.  This  is  the  reason  that  the 
prediction-error  power  of  .00837  above  is  not  zero.  When  we  increase  the  filter 
order  to  6,  the  results  of  table  B-2  and  figure  B-2  are  obtained.  We  notice  a 
significant  dc  indication  in  the  figure  and  a  nearly-zero  prediction  error  power 
of  -2.1E-10;  round  off  error  has  led  to  this  very  small  negative  number,  corresponding 
status  indicator  4,  and  a  negative  variance  from  the  graph.  We  also  notice  that  for 
this  noise-free  case,  all  six  poles  are  now  on  the  unit  circle  and  have  exactly 
the  correct  frequency,  except  for  round-off  error.  The  last  pole  is  actually  just 
outside  the  unit  circle  near  z=l,  and  corresponds  to  the  created  dc  component  noted 
earlier.  The  power  content  and  bandwidth  outputs  are  spurious  for  this  example, 
due  to  round-off  error. 

Lastly,  we  consider  the  effect  of  not  subtracting  off  the  sample  mean,  but 
keeping  p  at  5,  the  results  of  which  are  given  in  table  B-3  and  figure  B-3.  We 
again  see  perfect  pole  location  estimation  ,  as  predicted  by  (16) -(18),  and  a 
negative  prediction-error  power  of  -3.45E-11,  due  to  round-off  error.  The  spectral 
plots  in  figures  B-l  through  B-3  are  very  similar,  except  for  the  dc  indication  in 
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figure  B-2,  for  reasons  indicated  above. 

When  noise  is  added  to  the  tones,  the  significance  of  removing  or  leaving 
the  sample  mean  in  the  data  becomes  far  less.  It  appears  that  the  only  time  it  is 
important, is  when  the  true  dc  component  is  very  large,  in  which  case  the  sample 
mean  should  be  removed  so  as  not  to  contaminate  the  spectral  estimate  in  the 
neighborhood  of  zero  frequency.  All  the  results  in  the  main  body  of  this  memorandum 
have  had  the  sample  mean  removed. 
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LINEAR  PREDICTIVE  SPECTRAL  ANALYSIS  FOR  COMPLEX  DATA. 
FORWARD  AND  BACKWARD  AVERAGING.  MARPLE  ALGORITHM,  3  DEC  79. 
FILTER  WEIGHTS;  SPECTRUM;  POLE  LOCATIONS  AND  POWERS. 

NUMBER  OF  COMPLEX  DATA  POINTS  *  1QO 
MAXIMUM  ORDER  OF  FILTER  *  5 


MEAN  OF  DATA  =  -2. 00088008267E-02  -3. 7272203 195  IE-83 
VARIANCE  OF  DATA  =  4.37958610779 


STATUS  OF  STOPPING  CRITERION  =  1 
PREDICTION-ERROR  POWER  =  8. 37016559379E-03 
OPTIMUM  ORDER  OF  FILTER  =  5 


PREDICTION-FILTER  WEIGHTS: 


792285088326 
— « < 3633337 3 O  4 
-. 312567503922 
-. SO  1 325O04587 
-. S83323945939 


. 073S99SS3937 
-.262304443791 
. 717233365238 
.522347752791 
.726735507267 


VARIANCE  FROM  GRAPH  =  1.76044711955 


LOCATIONS  OF  POLES:  REAL, IMAG, MAG, ARG 

-.453612913509  -.890878008423 

Power  content  of  pole  =  .970737957165 
Amplitude  of  tone  =  .978052170292 


PI 


. 999714182602 
3  dB  Bandwidth 


1.32 


. 64991145304 
007937708E-04 


. 363261571624  -. 929531300933 

Power  content  of  pole  =  . 967105663763 
Amplitude  of  tone  =  .98026470355 


*3  Cj  jl;  C*  C*  <£+  *3  fT 

3  dB  Bandwidth 


99307O5775 

225032E-84 


.7703439763  . 637496046839 

Power  content  of  pole  =  .944648910737 
Amplitude  of  tone  =  . 979975939701 


. 999915522625  . 220052021064 

3  dB  Bandwidth  =  5. 37799640595E-05 


-. 53786963772  .383430154109 

Power  content  of'  pole  =  .  9994961 15623 
Amplitude  of  tone  =  .981047207863 


.999184665342  .69992545393 

3  dB  Bandwidth  *  5. S9937747441E-84 


-.390208035527  . 453332769386 

Power  content  of  pole  =  ,994133376334 
Amplitude  of  tone  =  .978736291359 


. 999012653532  . 85085702917 

3  dB  Bandwidth  =  .6. 2S5S42531 47E-84 


Table.  B-l.  Mean  SwUrvac^ed  }  J»  -  T 


Best  Available  Copy 
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LINEAR  PREDICTIVE  SPECTRAL  ANALYSIS  FOR  COMPLEX  DATA. 
FORWARD  AND  BACKWARD  AVERAGING.  MflRPLE  ALGORITHM,  5  DEC  79. 
FILTER  WEIGHTS;  SPECTRUM;  POLE  LOCATIONS  AND  POWERS. 

NUMBER  OF  COMPLEX  DATA  POINTS  *  100 
MAXIMUM  ORDER  OF  FILTER  =  6 


MEAN  OF  DATA  =  -2. 00080000267E-82  -3. 7272203 19A1E-A3 
VARIANCE  OF  DATA  =  4.37958610779 


STATUS  OF  STOPPING  CRITERION  *  4 
PREDICTION-ERROR  POWER  =-2. 13350977671E-18 
OPTIMUM  ORDER  OF  FILTER  =  6 
PREDICTION-FILTER  WEIGHTS: 


.295355531979 
.056631543087 
. 424649333142 
-. 238376931211 
032356553201 
. 634547103345 


. 0796434651 62 
-.343054271582 
. 931306306637 
-. 193517751776 
.204535335696 
-.723968642649 


VARIANCE  FROM  GRAPH 


=-43.4339553489 


LOCATIONS  OF  POLES:  REAL,  IMAG.,  MAG,  ARG/PI 

-.453990499742  -.391006524193  1 

Power  content,  of  pole  =  0  3  dB  Bandwidth  =  9 

Amplitude  of  tone  =  . 969947470875 


.363124552692  -.929776435333  1 

Power  content  of  pole  =  0  3  dB  Bandwidth  =  0 

Amplitude  of  tone  =  .938342333532 


.65 


-.379999999997 


.770513242771  .637423939744  .999999999995 

Power  content  of  pole  =  0  3  dB  Bandwidth  =  3 

Amplitude  of  tone  =  .980810321796 

-.537785252234  .389816994334  1 

Power  content  of  pole  »  8  3  dB  Bandwidth  =  8 

Amplitude  of  tone  =  . 981 1641 16999 

-.391806524135  .453990499756  1 

Power  content  of  pole  =  0  3  dB  Bandwidth  =  8 

Amplitude  of  tone  =  . 970330660887 


.219999999999 


.  69999999999 


. 34999999999 


1.80000801274  -3. 6436794 1244E-89  1.88830001269 

Power  content  of  pole  =  4 . 03994794443E-84  3  dB  Bandwidth 

Amplitude  of  tone  =  1. 51325361755E-8S 


-2. 7513S857492E-89 
=-3.37370491 132E-89 


la  tie  3-2,  -  6 
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LINEAR  PREDICTIVE  SPECTRAL  ANALYSIS  FOR  COMPLEX  BATA. 
FORWARD  AND  BACKWARD  AVERAGING.  MARPLE  ALGORITHM,  5  DEC 
FILTER  WEIGHTS;  SPECTRUM;  POLE  LOCATIONS  AND  POWERS. 
NUMBER  OF  COMPLEX  DATA  POINTS  =  130 
MAXIMUM  ORDER  OF  FILTER  *  5 


MEAN  OF  DATA  =  -2. 00800000267E-82  -3 . 7272233 1 95  IE-33  "<jfc  Sakracted 
VARIANCE  OF  DATA  *  4.37999999992 


STATUS  OF  STOPPING  CRITERION  =  4 
PREDICTION-ERROR  POWER  =-3 . 4471 5439332E- 1 1 
OPTIMUM  ORDER  OF  FILTER  *  5 


PREDICTION-FILTER  WEIGHTS: 


-.794144430743 
-.  73746294-2076 
-.312313565611 
-. 6316935446 
-.634547105931 


. 079648473335 
-.263405739893 
.717901319815 
.524333279376 
. 728963627434 


VARIANCE  FROM  GRAPH  =-l . 03533935768E-0S 


LOCATIONS  OF  POLES:  REAL, I MAG, MAG , ARG^P I 


-.453990499742  - . 39 1 386524 1 96  1 

Power  content  of  pole  *3  3  dB  Bandwidth  *  0 

Amplitude  of  tone  =  . 973 1 63678898 

. 363124552691  -. 929776485385  1 

Power  content  of  pole  =  3  3  dB  Bandwidth  -  0 

Amplitude  of  tone  =  .980342833582 

.778513242777  .63742398975  1 

Power  content  of  pole  =  8  3  d B  B an d w i d t h  =  8 

Amplitude  of  tone  *  .938313321793 


- .  65 


.379999999997 


.  219999999989 


-.587735252291  .889016994384  1  .7 

Power  content  of  pole  =  0  3  dB  Bandwidth  *  W 

Amplitude  of  tone  =  .931164116992 

-.39133652419  .453990499752  1  *8 

Power  content  of  pole  =  3  3  dB  Bandwidth  =  3 

Amplitude  of  tone  =  . 971817604263 

Talk  8-3.  5a*?l«-  Me**  N<t  )  r  =  * 
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APPENDIX  C.  DETECTABILITY  CRITERION 


Suppose  a  received  waveform  xflr)  is  characterized  as  follows : 


c  t),  H, 

n  (t),  H.  I 


(C-l) 


where  5 ft)  is  a  deterministic  signal  waveform  and  ft  (6  is  a  stationary  noise  process 
with  spectrum  let  us  process  x(B  linearly  according  to  correlation  operation 


■  =  dt  w  ft)  x  , 


(C-2) 


where  weighting  w(^  is  to  be  selected  to  maximize  the  deflection  of  output  random 
variable  2,  Letting  m  and  <r*  be  the  mean  and  variance  of  z,  we  define  deflection 

d  =  (c'3) 


(C-3) 


for  processor  (C-2) .  The  optimum  weighting  is  found  to  be 

W  -  W  ’ 

with  corresponding  maximum  value 

Now  suppose  that  signal  is  narrowband,  so  much  so  that  the  bandwidth 

is®!2  is  narrower  that  the  finest  detail  of  6r„(F)  at  the  signal  center  fre- 


(C-5) 


quency  fg.  Then  (C-5)  becomes 


<c  ■  m 


(C-6) 


where  E  is  the  received  signal  energy;  this  result  is  typical  in  white  noise 

backgrounds.  This  latter  quantity  is  a  "local”  measure  of  signal- to-noise  ratio 

at  f  ,  and  is  the  governing  statistic  for  detectability.  Values  of  (C-6) greater 
s 
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than  10  are  typically  needed  for  good  detection  and  false  alarm  probabilities. 
Let  us  now  consider  a  sampled  process,  with  time  increment  and 
(2a)  .  Then  we  can  let 


(C-7) 


where  J( u)  is  a  shape  function  defined  over  |u|<  I,  with  value  1  at  U  =  0.  Then 


noise  power 


%  I 

"  f  df  =■  du  J(u), 


(C-8) 


giving 


2*  tHjZt  f~ 


(C-9) 


in  terms  of  (T^  and  the  shape  function. 


If  our  signal  consists  of  a  pure  tone  of  frequency  f  and  amplitude 
A^,  and  duration  N  samples ,  then 

E  *  -y-A  M  >  (C-10) 


(C-10) 


and  (C-6)  yields 


A1  ,  AT/i  k, 

d“  c  Jtt/i) 


(C-ll) 


The  leading  term  in  (C-ll)  is  the  SNR  of  this  signal  tone,  referred  to  the  total 
noise  in  the  band  (-f^,f^) .  Hie  number  of  samples  available,  N,  enters  linearly 
into  this  result.  And  the  last  term  depends  on  the  shape  of  the  noise  spectrum. 

The  noise  example  we  will  consider  is  that  depicted  in  figure  14,  namely 


flv  =  3^  vv ,  +  .  4  w 


K-3  ,vk  > 


(C-12) 
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where  WK  is  white  noise  with  variance  6^ }  and 


=  or 7. 


(C-13) 


We  find,  for  example  (C-12) , 


'  - ^7* - /r  x  j^"'  > 

|  _  «;  e*f>(-ii J 


J  («)  = 


(\~  J 


C-  i  Trw) ) 


)  M<  I 


JduJto-  I±M±±g±S> 

Jo  (M-*) 


—  (  =  -+?. 2  dB  -p*r  *  =  o.-j) 


*.0P  CIS 


dbr  ®f  =  0  .y 


a _ t.  1  4-  4«?<2  +  °<* 

r”  ‘  *  R53 


(-Uis-ii  V  cf^t.n)  (c-\4) 


(C-15) 


(C-16) 


Equation  (C-15)  is  evaluated  in  table  C-l  for  the  parameters  of  figure 
22.  It  will  be  observed  that  values  of  d£  less  than  approximately  15  correspond 
to  tones  that  are  difficult  to  detect  in  figure  22,  whereas  values  of  d£  greater 
than  15  correspond  to  the  detected  tones  in  figure  22.  Thus  the  linear  predictive 
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technique  yields  not  only  the  spectral  estimate  of  the  process  under  investigation, 
but  also  detects  tones  at  local  SNRs  only  a  couple  of  dB  poorer  than  an  optimum 
processor.  The  quantitative  measure  in  (C-6)  is  a  very  useful  parameter  for 
quickly  ascertaining  the  detectability  of  tones  in  colored  noise. 


£ 

.07 

-10 

17.6 

.15 

-20 

13.2 

.35 

-30 

57.9 

.375 

-30 

80.6 

.45 

-30 

190.4 

.65 

-40 

89.1 

.75 

-50 

14.1 

.95 

-50 

22.0 

Table  C-l.  Deflections  for  Figure  22;  N  =  1000,  <  =  0.7 
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